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We have studied collisional perturbations of radio-frequency (rf) electric-dipole (El) transitions 
between the nearly degenerate opposite-parity levels in atomic dysprosium (Dy) in the presence of 
10 to 80 /iTorr of H2, N2, He, Ar, Ne, Kr, and Xe. Collisional broadening and shift of the resonance, 
as well as the attenuation of the signal amplitude are observed to be proportional to the foreign-gas 
density with the exception of H2 and Ne, for which no shifts were observed. Corresponding rates 
and cross sections are presented. In addition, rates and cross sections for O2 are extracted from 
measurements using air as foreign gas. The primary motivation for this study is the need for accurate 
determination of the shift rates, which are needed in a laboratory search for the temporal variation 
of the fine-structure constant [A. T. Nguyen, D. Budker, S. K. Lamoreaux, and J. R. Torgerson, 
Phys. Rev. A 69, 22105 (2004)]. 
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I. INTRODUCTION 

Perturbations of atomic resonance lines by collisions 
with atoms were first observed by Michelson 0. With 
the impact approximation, Lorentz 0] and Weisskopf Q 
provided the earliest theoretical framework. While 
Lorentz treated the effects of collisions as a complete ter- 
mination of a radiation wave train emitted by the atom, 
Weisskopf realized that a disruption in phase was suffi- 
cient to explain broadening. However, neither approach 
was able to explain shifts of the atomic resonances. Lind- 
holm J4| and Foley |5j remedied this shortcoming with the 
adiabatic impact approximation theory which takes into 
account the effects of distant collisions that induce small 
phase changes ignored by Weisskopf. While this the- 
ory is able to predict the line shape near the resonance, 
the conditions necessary for its validity break down near 
the wings or at higher pressures where satellite lines and 
asymmetries have been observed. The initial theory for 
satellite lines, called the quasistatic approach, was devel- 
oped by Khun |(J. In this approach, the motion of the 
perturber is ignored. The line shape is obtained by av- 
eraging the change in the energy difference between the 
atomic levels due to the presence of the perturber over 
the probability distribution of the stationary perturbers. 
However, this theory breaks down for frequencies near 
resonance or at low pressures. Since then, these failures 
have been alleviated by full quantum-mechanical treat- 
ments, developed to unify these two approaches. These 
techniques and related experiments are reviewed, for ex- 
ample, by Ch'en and Takeo 0, Allard and Kielkopf 
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and Szudy and Baylis 0. 

In this study, we investigate collisional perturbations 
of radio-frequency transitions between the nearly degen- 
erate opposite parity levels of 162 Dy and 163 Dy (Z=66) 
due to the introduction of H2, N2, He, Ar, Ne, Kr, and Xe 
at pressure ranging from 10 to 80 /iTorr. While broad- 
ening and shift rates, and corresponding cross sections 
are presented for 162 Dy, only shift results are presented 
for 163 Dy due to large statistical uncertainties resulting 
from the dilution of the level population by hyperfme 
levels. Rates and cross sections due to O2 are extracted 
from measurements using air as foreign gas. In addition, 
upper bounds on collisional-quenching cross sections are 
extracted from observations of the signal amplitude as a 
function of foreign-gas pressure. 

The system under consideration is unusual for a 
collisional study since the foreign-gas pressures are 
< 80 /iTorr. In this pressure regime, the time between 
collisions is longer than the total transit time of the Dy 
atoms across the apparatus, assuming a typical cross sec- 
tion of 10 -14 cm 2 . The only other study of collisional per- 
turbations in Dy that we are aware of was performed with 
a magnetic trap where the spin relaxation of the ground- 
state Zeeman sublevels due to He was presented [lOf . 

The primary motivation for this study is the need for 
accurate determination of the shift rates which are im- 
portant for a search for the temporal variation of the 
fine-structure constant, a [Tl| . Recently, evidence for 
variation in a on cosmological time scales was discov- 
ered in quasar absorption spectra [p3 . Il^ | , correspond- 
ing to a/a = (6.40 ± 1.35) x 10~ 16 /yr assuming a linear 
shift over 10 10 years [l3|. The current terrestrial evidence 
comes from an analysis of fission products of a natural 
reactor in Oklo (Gabon) that operated 1.8 x 10 9 years 
ago, corresponding to a/a — (— 2.3^°^) x 10~ 17 /yr [l4| . 
However, observational measurements of this kind are 
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FIG. 1: Relevant levels and transitions in atomic dysprosium. Population and detection scheme: Level B is populated in a 
three-step process. The atoms are excited to level f in the first two steps using 833- and 699-nm laser light (solid arrows). The 
third step is a spontaneous decay (labeled by a short-dashed arrow) from level f to level B with ~ 30 % branching ratio. The 
population is transferred to level A by the rf electric field (curved long-dashed lines). Atoms in level A decay to level c and 
then to the ground state. The fluorescence from the second decay (564 nm) is used for detection. Inset: Hyperfine structure of 
levels A and B for 163 Dy. Zero energy is chosen to coincide with the lowest hyperfine component. 



difficult to interpret due to assumptions and systematic 
uncertainties. In fact, both the astrophysical evidence 
and the Oklo analysis are in disagreement with differ- 
ent quasar absorption measurements |l5i llfij ] and earlier 
Oklo analyses [13, Ha] i respectively. These results have 
sparked interest in laboratory searches (see, for exam- 
ple, Refs. SH 0), which are easier to interpret 
but place weaker limits. The best laboratory limit of 
\6t/a\ < 1-2 x 10~ 15 /yr was obtained from a comparison 
of a Hg + optical clock to a Cs microwave clock [2(j ■ 

Recent calculations [23 have shown that the nearly 
degenerate opposite-parity levels in atomic dysprosium 
(Fig. ^ levels A and B) are highly sensitive to varia- 
tions in a. In Ref. a possible ultimate experimental 
sensitivity of \a/a\ ~ 10" 18 /yr was estimated for this 
system from an analysis of statistical and systematic un- 
certainties. In order to realize this level of sensitivity, the 
transition frequencies must be kept stable to a mHz level 
over a period of a year. Thus, knowledge of collisional 
shift rates is critical for the determination of the maxi- 
mum allowable pressure in the vacuum chamber and, if 
necessary, for correction of the frequency measurements. 

In Section II, we describe the technique used to mea- 
sure transition frequencies and the procedure used dur- 
ing collisional-perturbation measurements. Section III 
presents the results and the details of the analysis. Sec- 
tion IV considers various systematic uncertainties. In 
Section V, we extract broadening and shift cross sec- 
tions. In Section VI we present the implications for the 
a- variation experiment. 



II. EXPERIMENTAL TECHNIQUE 

A. Overview 

The energy difference between the nearly degenerate 
levels considered in Refs. |2^| and ^lj are on the or- 
der of hyperfine splittings (Fig. ^ inset) and isotope 
shifts. Since Dy has seven stable isotopes, there are many 
rf transitions. The 235-MHz 162 Dy (4f 9 5d 2 6s J=10 -> 
4f 10 5d6s J=10) transition and the 3.1-MHz (F=10.5 -> 
F=10.5) transition are used for the a-variation experi- 
ment due to practical considerations such as high count- 
ing rate and low frequency, respectively. As shown in 
Fig. ^ the population transfer from the ground state to 
the long-lived (tb > 200 /xs |24|) odd-parity level B re- 
quires three transitions. The first two transitions are in- 
duced via 833- and 669- nm laser light, while the last tran- 
sition is a spontaneous decay with a ~ 30 % branching 
ratio [2{| . Atoms are transferred to the even-parity level 
A (ta=7.9 /its [24|]) with a frequency-modulated rf electric 
field referenced to a commercial Cs frequency standard. 
The transition frequency is determined via a lock-in tech- 
nique by monitoring the 564-nm fluorescence. 



B. Apparatus 

In this section, we describe the atomic-beam appara- 
tus, the optical setup, and the detection and pressure- 
measurement systems. 

A detailed description of the atomic-beam source is 




FIG. 2: Experimental setup (not shown to scale): a) Atomic 
beam produced by effusive oven source at ~ 1500 K; b) 
atomic-beam collimator; c) oven light collimator; d) cylin- 
drical lenses to diverge laser beams; e) interaction region of 
atoms with the electric field enclosed in a magnetic shield 
(not shown); f) cylindrical mirror to collect fluorescent light; 
g) lucite light pipe; h) interference filter; and i) short-pass 
filter. 



given in Ref. |25j . The beam is produced by an effusive 
oven with a multislit nozzle-array operating at ~ 1500 K. 
The oven consists of a molybdenum tube containing dys- 
prosium metal, surrounded by resistive heaters made 
from tantalum wire inside alumina ceramic tubes. The 
radiation shielding is provided by five layers of tantalum 
sheets surrounding the oven and the heaters. In addition 
to the multislit nozzle, two external collimators are used 
to collimate the atomic beam and the oven light. The lat- 
ter is necessary to minimize the background due to scat- 
tered oven light. The resulting atomic beam has a mean 
velocity of ~ 5 x 10 4 cm/s with a full-angle divergence of 
~ 0.2 rad (1/e 2 level) in both transverse directions. 

The atoms enter the interaction region (Fig. [21, where 
they are first transferred to level B by two laser-induced 
transitions and a spontaneous decay. Approximately 
600 mW of 833-nm light is produced by a Ti:Sapphire 
ring laser (Coherent 899) pumped by 15 W of Ar ion laser 
light (Coherent Innova 400). A ring dye laser (Coherent 
699) with DCM [4-(Dicyanomethylene)-2-methyl-6-(p- 
dimcthlaminostyryl)-4H-pyran] dye produces ~ 300 mW 
of 669- nm light. The dye laser is also pumped by an 
Ar ion laser (Coherent Innova 300) operating at 6 W of 
power. 

Since the population transfer of a weakly collimated 
atomic beam due to narrow-band cw lasers is inefficient, 
an adiabatic passage technique is utilized to transfer 
atoms to level B. A detailed description of this tech- 
nique and references to earlier work are given in Ref. |26j . 
Briefly, cylindrical lenses are used to diverge the two laser 
beams in the interaction region such that the divergences 
of the light beams match the atomic-beam divergence. 
Due to the Doppler effect, the atoms "see" a frequency 
chirp in laser detuning which adiabatically transfers the 
population to the excited state. 
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FIG. 3: Schematic diagram of the rf-generation, control, and 
pressure-measurement systems. The rf synthesizer is refer- 
enced to a Cs standard, which is compared to a CDMA (Code 
Division Multiple Access) primary reference source. Both 
clocks have a long term stability of ~ lO -11 . The output of 
the synthesizer is amplified and then sent to the electric-field 
plates and a 50 Q power-resistor termination in series with 
the larger impedance of the plates. The output of the pho- 
tomultiplier tube is sent to the lock-in amplifier for detection 
and acquisition. 



The rf electric field is formed between two parallel, 
10 cm by 5 cm, wire grids made of 0.002 in Be-Cu 
wire. The separation between adjacent wires is 2.5 
mm which minimizes surface area and hence surface- 
charge accumulation. The separation between the grids is 
~ 2.5 cm. The entire assembly is enclosed inside a high- 
permeability magnetic shield which reduces the magnetic 
field in the interaction region to ~ 1 mG. The 564-nm flu- 
orescence is collected by a lucite light pipe and detected 
with a photomultiplier tube (PMT) with bandpass inter- 
ference and short-pass filters at the input window. 

The rf-generation and detection system is shown in 
Fig. [3J The frequency-modulated rf field is generated by 
a synthesizer (Hewlett Packard 8647A) and sent through 
a 45-dB amplifier (ENI 604L). The rf power at the out- 
put is ~ 0.5 W. The synthesizer is referenced to a Cs 
frequency standard (Hewlett Packard 5061A) with a long 
term stability (>100 s) of 10 -11 . This frequency standard 
is compared to a second clock (Symmetricom TS2700) 
to monitor its stability. A Signal Recovery 7265 lock-in 
amplifier provides a 10-kHz modulation signal to the rf 
synthesizer and demodulates the signal from the PMT. 
The signal from the lock-in along with the clock, pres- 
sure, and temperature data are sent to a computer for 
data acquisition. 

The pressure measurement is provided by a residual 
gas analyzer (SRS RGA 200) and a capacitance manome- 
ter (MKS Baratron 690A.01). The residual gas analyzer 
(RGA) is used to monitor partial pressures of the back- 
ground gasses. However, due to unknown gas-species- 
dependent sensitivity factors, it is not well suited to 
measuring the absolute pressure of various gases. The 
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manometer, on the other hand, is capable of absolute 
pressure measurement in a wide range of pressures from 
0.1 Torr to a few /iTorr with a resolution of 10~ 7 Torr. 
Measurements have shown that the manometer response 
is nonlinear below 10 /xTorr. Thus, the total manometer 
pressure was limited to 20 /xTorr during the collisional- 
perturbation measurements. The manometer signal is 
processed by a signal conditioner (MKS 670B) which 
includes an active heater control mechanism to stabi- 
lize the operating temperature. To further reduce noise, 
the manometer is mechanically isolated by a vibration- 
isolation mount at the base and by the use of flexible 
bellows to connect to the vacuum chamber. The mea- 
sured statistical uncertainty in the pressure range used 
during measurements is ~ 1 /iTorr for 100 seconds of 
integration. 



C. Procedure 

To minimize fluctuations such as those due to laser- 
power drifts or density fluctuations of the atomic beam, 
the rf voltage is frequency modulated at 10 kHz with an 
equal modulation depth (modulation index of 1). Since 
the rf-transition linewidth is ~ 20 kHz, this modulation 
provides a fast sweep across the absorption line shape. 
The outputs from the lock-in amplifier at the first and 
second harmonic of the modulation frequency are shown 
on Fig. 0] The first-harmonic signal is an odd function 
of detuning while the second-harmonic signal is an even 
function. It is well known that for small modulation am- 
plitudes, the output signal at the harmonics of the mod- 
ulation frequency are proportional to the corresponding 
derivatives of the absorption line shape j^. In our sys- 
tem, since the modulation amplitude is comparable to 
the rf-transition linewidth, the harmonic signals contain 
higher derivative contributions. Therefore, a fitting func- 
tion was derived from first-order time-dependent pertur- 
bation theory. This model does not include any effects 
due to power broadening or beam velocity distribution 
and is used to confirm qualitative behavior of the line 
shape. 

To measure the transition frequency, the ratio of the 
first- and second-harmonic signals is used in order to 
reduce drifts further. The transition frequency is ex- 
tracted from the ratio by a two-step process. First, 
the lock-in amplitude at the first and second harmon- 
ics and their ratio are recorded simultaneously |28| as 
the carrier frequency is scanned 4 kHz around the reso- 
nance. For such a small scan compared to the linewidth, 
the first harmonic is linear while the second harmonic 
is nearly a constant (Fig. 0}. By fitting a linear func- 
tion to the ratio, a conversion parameter from the ratio 
to frequency is obtained. Thus, these scans are called 
calibration files. Once the calibration is established, the 
rf carrier frequency is locked at a fixed value near res- 
onance (3,074,000 Hz and 234,661,000 Hz for the two 
transitions) , and the ratio is measured repeatedly. These 
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FIG. 4: First- and second-harmonic signals for 10-kHz fre- 
quency modulation. Solid line is a fit to the function derived 
from first-order time-dependent perturbation theory. 

measurements are saved in the so-called fixed-frequency 
files. The average of ratios from these measurements are 
converted to frequency using the calibration file. 

It should be mentioned that before recording a calibra- 
tion file, a specific phasing procedure is used to maximize 
the lock-in-amplifier signal at both harmonics. Since 
Signal Recovery 7265 lock-in amplifier allows for sepa- 
rate phase adjustments for the two harmonics, the first- 
harmonic signal phase is adjusted to maximize the sig- 
nal at 10 kHz above the fixed-frequency value, while the 
second-harmonic signal is maximized at the fixed fre- 
quency value. This procedure is used since the first- and 
second-harmonic signals are maximal near these frequen- 
cies, respectively. In addition, to minimize any instability 
or loss of signal due to laser detuning, each laser is re- 
tuned between recording data files by scanning over the 
laser resonances and setting the laser frequency to the 
resonance peak. Using this technique for frequency mea- 
surements, we have been able to achieve ~ 2- and 6-Hz 
statistical sensitivity for 1 second of integration time for 
the 235- and 3.1-MHz transitions, respectively. 

For collisional-perturbation measurements, foreign gas 
is introduced into the chamber through a leak valve (Var- 
ian Model 1000). In order to minimize contamination, 
the hose connecting the gas cylinder to the leak valve 
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FIG. 5: Shift of the resonance frequency for both transitions 
due to the presence of Krypton. Av is the difference between 
the transition frequency and the fixed rf synthesizer carrier 
frequency (3,074,000 Hz and 234,661,000 Hz for the two tran- 
sitions) . The solid lines are the least-squares linear fits to the 
data. 

is evacuated with a separate pump and flushed with gas 
several times prior to opening the valve. Then, the leak 
valve is opened and foreign gas is introduced continuously 
until a dynamic equilibrium is reached at a manometer 
pressure reading of ~ 80 /iTorr. Meanwhile, the RGA is 
monitored for contamination due to other residual gases 
in the hose. After the background gas pressures drop 
below 10 /iTorr, measurements are taken at descending 
pressures between 80 and 20 /iTorr. At each pressure, one 
calibration and three fixed-frequency files are taken for 
each transition. To check for systematic drifts and signal 
loss, two more measurements at intermediate pressures 
are taken after the lowest-pressure measurement. 

Since foreign gas is introduced from the top of the 
chamber while the chamber is continuously pumped from 
the bottom, a pressure gradient between the top and the 
interaction region might exist. Due to the large conduc- 
tance of our chamber, we estimate a pressure gradient 
of ~ 5%. This estimate is in agreement with the upper 
bound extracted from measurements taken with varying 
pumping speeds. Since the uncertainty in this pressure 
calibration error is not well known, this correction is not 
included in the results presented in Sec 11111 and Sec. El in 
order to preserve the relative precision of the measure- 
ments. 



III. ANALYSIS AND RESULTS 

A. Frequency shifts 

For all the gases except H2 and Ne, which do not shift 
the transition frequency, ~ 100 Hz shifts were observed 
over the 10 to 80 /iTorr pressure range. Since these tran- 
sitions are between two excited states, the measured shift 
rates are due to the difference between the shift of the 
levels involved in the transition. As shown in Fig. \5\ 



TABLE I: Collisional-shift rates for the 235-MHz and 3.1- 
MHz transitions in the presence of various gases. The shift 
rates are expected to be equal and opposite for the two tran- 
sitions since there is a change in sign between the energy 
difference of the levels involved. The O2 coefficient was ex- 
tracted using a measurement when air was let into the vacuum 
chamber. This was only performed for the 235-MHz transi- 
tion. The results have been corrected for the presence of the 
background gases and the assigned errors include the effects 
of Dy self-shifts discussed in Sec. If VI As discussed at the end 
of Sec. Ill CI the analysis does not include the estimated 5% 
correction due to pressure calibration error. 

235-MHz Transition 3.1-MHz Transition 

shift rates (Hz//iTorr) shift rates (Hz//iTorr) 



H 2 -0.02 ± 0.02 -0.09 ± 0.06 

N 2 -1.71 ± 0.03 1.72 ± 0.05 

2 -1.97 ± 0.25 

He 1.25 ± 0.02 -1.27 ± 0.05 

Ne -0.01 ± 0.02 -0.02 ± 0.05 

Ar -2.21 ± 0.05 2.14 ± 0.07 

Kr -2.78 ± 0.05 2.78 ± 0.07 

Xe -2.74 ± 0.04 2.75 ± 0.07 



the shifts for the two rf transitions are in opposite direc- 
tions. This is expected since for the 235-MHz transition, 
level B is lower in energy than level A, while for the 3.1- 
MHz transition, the opposite is the case. Moreover, since 
electronic structures of the isotopes are identical, the ab- 
solute value of the shift coefficients should be identical, 
which is indeed the case within the precision of the ex- 
periment. For the gases that produce shifts, all but He, 
shift a given transition in the same direction. 

The analysis of the shift data is complicated by the 
effects of background gases. Since the normal operating 
pressure of the vacuum chamber with the oven turned 
on is ~ 5 - 10 /iTorr, the background gases, dominated 
mostly by H2 and N2, and to a lesser extent by O2, make 
up a significant fraction of the overall pressure. More- 
over, in some cases these gas pressures are not constant 
throughout a measurement |29| and contribute to the ob- 
served shift rates. 

Since the capacitance manometer measures the total 
gas pressure in the chamber, RGA partial gas pressures 
are needed to establish the dominant-gas pressure and 
correct for the systematic shifts induced by the back- 
ground gases. However, RGA partial pressures must be 
corrected with gas-species-dependent sensitivity factors 
that are not well determined. We have tried to extract 
these sensitivity factors using both data taken during 
collisional-perturbation measurements, where the oven is 
operating at ~ 1500 K, and auxiliary cold-oven measure- 
ments, where gas is introduced into the chamber with 
the oven turned off. In the latter case, the large H2 
background is absent, and the operating pressure of the 
vacuum chamber (without the foreign gas) is ~ 1 /iTorr, 
dominated by N 2 . Therefore, the background gases can 
be ignored and the sensitivity factors extracted using lin- 
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ear fits to the manometer versus RGA data. 

Extracting these parameters from the hot-oven data 
relies on the fact that during measurements where H2 is 
introduced into the chamber, the N2 and O2 pressures 
are stable (we ignore all other gases including H2O since 
their pressures are below 5 x 10~ 7 Torr throughout these 
measurements). Using this fact, the H2 sensitivity fac- 
tor, Ph2, was extracted in the same way as the auxil- 
iary cold-oven method, ignoring the other gases. The 
N 2 sensitivity factor, 13^2, was extracted from measure- 
ments where N 2 was introduced into the chamber, but 
by accounting for the background H2 pressure (using the 
extracted (3h2)- Finally, the O2 sensitivity factor, 0O2, 
was extracted taking both the N2 and H2 pressures into 
account using data where air was introduced into the 
chamber. 

The RGA sensitivity factors extracted using the cold- 
oven and hot-oven data disagree at ~ 10 - 20 % level, sug- 
gesting that they are also dependent on the background- 
gas pressures. As a result, in the rest of the analysis, 
sensitivity factors extracted from in-situ hot-oven mea- 
surements were used with larger error bars to accommo- 
date for the results obtained by the auxiliary cold-oven 
measurements. In addition, in order to minimize the un- 
certainty introduced by the use of RGA sensitivity fac- 
tors, the dominant gas pressures, P gas , were calculated 
using 

P gas = P(CM) - faPi(RGA), (1) 

i 

where P(CM) is the manometer pressure, Pi(RGA) and 
Pi are the RGA partial pressure and sensitivity factor 
of a background gas species i, respectively. Since the 
background-gas pressures are at least two times smaller 
than the dominant-gas pressure, by using this method, 
the uncertainty due to the sensitivity factors are reduced 
to levels comparable to or smaller than the manometer 
uncertainty of 1 /iTorr. 

Once the sensitivity factors of the background gases 
were determined, an iterative approach was used to ex- 
tract the Dy rf-transition-shift rates for H2, N2, and 02- 
The iterative approach is necessary for these gases be- 
cause for a measurement where one of these gases domi- 
nates, the other two are background gases. As a result, 
to extract the correct shift rate for one of them requires 
the knowledge of the shift rates of the other two, which, 
in turn, are dependent on the shift rate that is being ex- 
tracted. So, on the first iteration, the H2 shift rate was 
extracted by taking into account the background N2 and 
O2 pressures but not the shifts due to them. Then, the N2 
shift rate was extracted similarly, by taking into account 
both the H2 and O2 pressures, as well as shifts due to H2. 
Finally, the O2 shift rate was extracted by taking into ac- 
count both N 2 and H 2 pressures and shifts (only for the 
235-MHz transition since no air data were taken for the 
3.1-MHz transition). Then, these first-order coefficients 
were used to extract higher-order coefficients where pres- 
sures and shifts due to all the background gases were used 
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FIG. 6: Collisional broadening for the 235-MHz transition 
due to Kr without the phasing-effect correction. On the y 
axis, the inverse of the normalized first harmonic slope, which 
is directly proportional to the linewidth is plotted. The solid 
line is the least-squares linear fit to the data. 

in every step. This iterative process converged within 
three iterations, improving x 2 at every step. 

The rest of the gases were easier to analyze since they 
do not appear as background gases in other measure- 
ments. The coefficients for these gases were determined 
by taking into account the background H2, N2 and O2. 
For the 3.1-MHz transition, the O2 shift rate for the 235- 
MHz transition with the opposite sign was used. The 
results are presented in Tabled 



B. Broadening 

Since the main focus of this study is frequency shifts, 
the experimental technique we utilize is optimized for 
frequency measurements. Nevertheless, the broadening 
rates can also be extracted from the calibration files. 
However, determination of the broadening rates is more 
difficult than the shift rates since, as discussed below, 
the sensitivity of the calibration data to linewidth is low. 
In addition, there is a systematic effect caused by the 
phasing procedure used (see Sec. Ill G|) coupled with the 
collisional shift of the resonance frequency that results 
in an enhancement of apparent broadening rates for one 
transition and reduction for the other. 

The slope of the first-harmonic signal near resonance, 
normalized by the second harmonic can be extracted by 
fitting a linear function to the ratio from the calibration 
file data, as described in Sec. Ill Gl For small modula- 
tion amplitudes, the first-harmonic signal is proportional 
to the first derivative of the absorption line shape and 
the second harmonic to the second derivative. Near reso- 
nance, to first order in frequency detuning, both of these 
functions have a I/7 3 dependence. As a result, the slope 
of the ratio is insensitive to the linewidth. However, in 
our system the modulation amplitude is comparable to 
the linewidth, and the harmonics contain higher-order 
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TABLE II: Collisional broadening rates for the 235-MHz 
transition in the presence of various gases. The large uncer- 
tainties are due to statistical noise and the uncertainty in the 
conversion factor from the inverse slope to linewidth. The O2 
result is presented as an upper bound at 2a confidence level. 
The results do not include the estimated 5% uncertainty due 
to pressure gradients in the chamber. 

235-MHz Transition 
broadening rates (Hz//iTorr) 



H 2 


8.6 ± 1.9 


N 2 


3.0 ± 1.7 


2 


< 8.6 


He 


5.3 ± 1.9 


Ne 


1.4 ± 1.3 


Ar 


4.9 ± 1.4 


Kr 


5.3 ± 1.4 


Xe 


2.8 ± 2.0 




4x10 
Pressure (Torr) 

FIG. 7: Decrease in the second harmonic signal for the 235- 
MHz transition due to Xe . The signal has been normalized 
to one at zero pressure. The solid line is the least-squares fit 
to the data. 



derivatives which result in a small dependence of the 
slope on the linewidth. In practice, the inverse of this 
slope is the useful quantity in converting the ratio to 
frequency. It is proportional to the linewidth for small 
changes in the linewidth 30] and increases with gas pres- 
sure due pressure broadening (Fig. . 

The conversion factor from the inverse slope to 
linewidth is extracted from measurements where calibra- 
tion files and line-shape scans were taken simultaneously 
for different rf input powers, which lead to power broad- 
ening. The extracted conversion factor of 31(5) agrees 
with the conversion factor of 28.2 obtained from the first- 
order perturbation theory calculation by linearizing the 
ratio near resonance. In the analysis of the data, the ex- 
perimental value is used and errors propagated accord- 
ingly. 

There is a systematic effect related to the phasing 
procedure which is due to the dependence of the first- 
harmonic slope near resonance, and hence the slope of 



the ratio, on the frequency at which the phase is adjusted 
to maximize the signal (which we define as phasing fre- 
quency) . For small changes in the phasing frequency, the 
experimental measurements show a linear dependence 
with 0.025(3) and 0.029(5) Hz change in the inverse slope 
per Hz change in the phasing frequency for the 235- and 
3.1-MHz transitions, respectively. For all the gases where 
shifts are observed, with the exception of He, which shifts 
the transition frequency in the opposite direction, this 
effect leads to an enhancement of the broadening rates 
for the 235-MHz transition and a reduction for the 3.1- 
MHz transition. This is due to the fact that as the gas 
pressure is changed, the transition frequency shifts but 
the phasing frequency is fixed at 10 kHz above the carrier 
frequency used for the fixed-frequency measurements. As 
a result, the phasing frequency changes with respect to 
the transition frequency. For example, from the Kr shift 
data shown in Fig. El at low pressures, the rf carrier fre- 
quency for the 235-MHz transition is 50 Hz below the 
transition frequency, while at high pressures, it is 100 Hz 
above the transition frequency (the y axis of the figure is 
transition frequency minus carrier frequency). As a re- 
sult, while at low pressures, the phasing frequency with 
respect to the transition frequency is less than 10 kHz, 
at high pressures, it is greater than 10 kHz. This leads 
to an enhancement of apparent broadening since at low 
pressures the observed inverse slope is less than the value 
that would be observed without this effect, and the oppo- 
site is the case at high pressures. The 3.1-MHz transition 
is shifted in the opposite direction, and as a result, the 
observed broadening is reduced for this transition. This 
is a significant effect since the changes in the transition 
frequency lead to 3 - 4 Hz changes in the inverse slope. In 
comparison, as seen in Fig. [fJl the changes in the inverse 
slope due to the gases are 10 - 15 Hz. 

The results that have been corrected for this systematic 
effect, as well as the systematic effect due to the presence 
of the background gases, arc presented in Table ITT1 Even 
with this correction, the statistical noise coupled with 
the uncertainty in the conversion from inverse slope to 
linewidth lead to poor fits for the 3.1-MHz data. Conse- 
quently, the broadening results for this transition are not 
presented. For the same reasons, only an upper bound 
for the 235-MHz O2 broadening rate is presented. Due to 
small O2 pressure variation (< 1.1 /iTorr), the omission 
of the O2 correction is not problematic in this case since 
even if the broadening rate for O2 were 8.6 Hz//iTorr 
(the upper 2a bound from Table IXTf) . this would lead to 
an error of < 5% for the measurements that exhibit the 
largest O2 pressure variation. 



C. Signal Attenuation 

During the fixed-frequency measurements, ~ 15 - 40 % 
decrease in the second-harmonic-signal amplitude, which 
is nonzero at resonance, as a function of foreign-gas pres- 
sure was observed (Fig. 0). This signal attenuation can 
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TABLE III: Signal-attenuation rates for the 235-MHz and 3.1-MHz transitions and upper bounds for collisional quenching 
cross sections in the presence of various gases. The results have been corrected for the presence of the background gases 
and the assigned error bars include systematics due to beam-density variations. The results do not include the estimated 5% 
uncertainty due to pressure gradients in the chamber. For the combined cross sections, gas temperature of 300(5) K and Dy 
beam temperature of 1500(50) K are assumed. 





235-MHz Transition 
e (10 -3 1/^iTorr) 


3.1-MHz Transition 
e (10^ 3 l//iTorr) 


Quenching 

<a q v> (10~ 10 cm 3 /s) 


H 2 


4.1 ± 0.2 


3.9 ±0.2 


< 13 


N 2 


4.4 ±0.5 


4.3 ±0.7 


< 14 


o. 


10 ± 1 




< 32 


He 


2.6 ± 0.2 


2.6 ± 0.2 


< 8 


Ne 


2.3 ± 0.2 


2.6 ± 0.2 


< 7.6 


Ar 


3.3 ±0.3 


3.7 ±0.2 


< 11 


Kr 


3.9 ±0.2 


4.0 ± 0.2 


< 13 


Xe 


4.9 ± 0.2 


5.2 ± 0.2 


< 16 



be attributed to either decrease in the beam density due 
to elastic scattering or decrease in the fluorescence in- 
tensity due to inelastic quenching. The signal amplitude 
decreases with the distance the beam travels through the 
foreign gas: 



<S(v)>= GI a <e 



-l/\(v) 



>, 



(2) 



where G is the efficiency of the detection system, Iq is 
the fluorescence intensity at zero pressure, A is the mean 
free path, and / is the beam path length through the 
foreign gas. Here, < ... > denotes average over the beam 
velocity distribution. Since A is inversely proportional to 
pressure, the signal for small pressures can be written as: 



<S{v)> 
So 



=< e -«A(*) >RJ <i _ i/ X (v) >= (l - € p), (3) 



where So is the second-harmonic signal at zero pressure, 
and e is the attenuation rate. 

The analysis of the quenching data is similar to the 
analysis of the shift and broadening. First, the iterative 
process was used for the dominant gases in the back- 
ground, H2, N2, and 02- Then, the values extracted 
were used in the analysis of the remaining gases. The at- 
tenuation rates tabulated in Table lull include systematic 
corrections due to clogging of the unheated photon colli- 
mator (Fig. and temperature drift of the ovens. To de- 
tect any changes in the signal amplitude due to clogging, 
the two measurements at intermediate pressures taken 
after the lowest pressure points, as discussed in Sec. Ill CI 
were used. The oven temperature was monitored by a 
thermocouple and changes in the beam density due to 
temperature fluctuations were estimated using Dy sat- 
urated vapor-pressure dependence on temperature (see 
Sec. US) 

The interpretation of the signal attenuation data is 
complicated since it is difficult to separate elastic scat- 
tering events that decrease beam density from inelastic 
quenching that decreases the population of level B. Nev- 
ertheless, we expect the elastic scattering contribution to 



be small since although typical elastic cross sections are 
~ 10~ 14 cm 2 (corresponding to a mean free path of ~ 30 
cm at 100 /iTorr), large angle scattering necessary for sig- 
nificant collimation loss is suppressed by more than one 
order of magnitude due to the fact that differential cross 
section is sharply peaked at small scattering angles [3l| . 
Using the results derived in Ref. f° r hard-ball scat- 
tering, we have estimated that the differential cross sec- 
tion becomes appreciable below 0.1 - 0.01 radian in the 
center-of-mass frame. For a head-on collision, this results 
in 6 - 7 mrad deviation for the Dy beam in the laboratory 
frame. This deviation is small compared to the collima- 
tion half angle of 0.1 radian used in our experiment. How- 
ever, since this estimate uses a hard-ball model, which 
depends on few assumptions such as a specific minimum 
distance between the centers of the spheres, as opposed 
to more complicated interactions such as quasi-molecular 
formation, it can vary by a multiplicative factor of 2 - 5. 
As a result, we only present upper bounds for the inelas- 
tic quenching cross sections. 

The attenuation rates can be related to the quenching 
cross sections through a characteristic mean free path, 
which can be expressed as 



A 



<VB> 
no <<JqV> 



(4) 



where <v B >= 3/4(2?r k T B /m B ) 1/2 is the average beam 
velocity given by a modified Maxwellian distribution, no 
is the density of foreign gas, v is the relative velocity 
of the gas and beam particles, and <a q v> denotes an 
average over both the gas and the beam velocity distri- 
butions. 

Using Eq. (0} and Eq. the quantity <a q v> can be 
written in terms of e as 



<<J n V>= -t 



-TgT 1 / 2 , 

TODy I 



(5) 



where Tq and T are the temperatures of the foreign gas 
and the oven, respectively. The path length, I, is not well 
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defined in our experimental setup. Since the light pipe 
has a width of 6.7 cm, it collects fluorescence from atoms 
at various path lengths. The tabulated results assume a 
path length of 5.5 cm, which is the distance between the 
laser excitation region and the center of the light pipe. 
The results are summarized in Table ITTT1 



IV. SYSTEMATIC UNCERTAINTIES 

In addition to systematic effects considered above, 
Dy*(excited Dy atoms)-Dy* and Dy*-Dy collisions 
within the beam can lead to a systematic change in the 
transition frequency and linewidth, if the beam density 
does not remain constant during the collisional measure- 
ments. We have performed independent measurements of 
self-collisions in our system and present here an estimate 
of their contribution to the systematic uncertainty. 

The beam density can be controlled by varying the 
oven temperature. The temperature dependence of Dy 
saturated- va por pressure is given by the following empir- 
ical formula |33| : 



p _ 2Ql2.460-15336/T-1.1141og(T) 



(6) 



where P is the pressure in Torr and T is the temperature 
in Kelvin. 

We have varied the oven temperature, which was 
monitored by a thermocouple, from 1077°C to 1170°C, 
while monitoring the frequency for both transitions. 
This temperature range corresponds to a Dy density of 
(0.3 — 1.5) x 10 14 cm -3 in the oven. Earlier absorption 
measurements |25j have shown that the density of atoms 
in the population and detection region is ~ 5 x 10 9 cm -3 
at the normal operating temperature of 1170°C. This 
density is in agreement with estimates based on effusive 
flow. To obtain the Dy density in the interaction re- 
gion for various temperatures, the scaling factor between 
the calculated density in the oven at 1170°C and the ab- 
sorption measurements was applied uniformly to all the 
calculated densities. The extracted self-shift rates were 
0.2(2.0) x 10~ 10 Hz cm 3 and -3.1(6.0) x 10~ 10 Hz cm 3 
for the 235-MHz and 3.1-MHz transitions, respectively 
(Fig. (SJ . An upper bound for the magnitude of the self- 
shift rate can be estimated by noting that the rates, m, 
for these two transitions should be the same in magnitude 
but opposite in sign: 



M = ^l TO 235 ~ m 3 .i|, 
= |m 2 35 + w 3 .i|, 



(7) 



where the second equation provides a way to estimate 
errors. The resulting upper bound on the magnitude of 
the self-shift rate is < 4.9 x 10~ 10 Hz cm 3 at la confi- 
dence level. Since changing the oven temperature does 
not discriminate between Dy* and Dy densities, this up- 
per bound is an effective rate due to both Dy*-Dy* and 



-105 
-110 



2x10 4x10 
Dy density (cm -3 ) 



105 
100 




2x10 4X10 
Dy density (cm -3 ) 



FIG. 8: Shifts due to collision of the Dy atoms in the beam. 
Beam density was varied by changing the oven temperature 
from 1170°C to 1077° C which correspond to the estimated 
atomic densities in the interaction region shown in the plots. 



Dy*-Dy collisions. A simple way to address this ambigu- 
ity would be to change the laser powers simultaneously 
to change the fraction of excited Dy atoms. This method 
may be pursued in a future study. 

Using the extracted upper bound, we can estimate the 
systematic contribution of self-collisions to the measured 
shift and broadening rates. This requires a mechanism 
that changes Dy beam density as a function of foreign-gas 
pressure. As a worst-case scenario, we assume that all of 
the observed signal attenuation is due to the decrease in 
the Dy density as a result of scattering with the foreign- 
gas atoms. Depending on the gas, this results in a change 
of ~ (0.9 - 2.2) x 10 9 cm" 3 in the Dy density, which 
corresponds to ~ 0.4—1.1 Hz change in the magnitude of 
the transition frequency from low to high has pressures. 
The correction to the foreign-gas shift rates are in the 
range of 0.005 — 0.01 Hz//iTorr. These uncertainties have 
been included in Table [I] 

Since, in general, broadening should be of the same 
order of magnitude as shifts, the uncertainty due to self- 
broadening is negligible compared to the errors quoted 
on Table Oil 
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TABLE IV: Collisional shift and broadening cross sections and their ratios at foreign-gas temperature of 300(5) K. The shift 
cross sections for the two transitions have been combined into a single result with the following sign convention: a positive 
frequency shift for the 3.1-MHz transition gives a positive cross section. 

<a s v> (10~ 10 cm 3 /s) <a b v> (10 -10 cm 3 /s) <a s v>/2 <a b v> 



H 2 0.05 ± 0.04 8.4 ± 1.9 0.003 ± 0.002 

N 2 3.34 ±0.07 2.9 ±1.7 0.6 ± 0.3 

2 3.9 ±0.5 < 8.4 < 0.2 

He -2.44 ± 0.05 5.2 ± 1.9 -0.2 ± 0.1 

Ne 0.02 ±0.04 1.4 ±1.3 0.01 ± 0.02 

Ar 4.3 ±0.1 4.8 ±1.4 0.4 ±0.1 

Kr 5.4 ±0.1 5.2 ±1.4 0.5 ± 0.1 

Xe 5.3 ±0.1 2.7 ±2.0 1.0 ± 0.7 



V. DISCUSSION 

The adiabatic impact approximation results give a 
Lorentzian line shape with the transition frequency and 
linewidth modified by |34|: 



. 1 , n . <va b > 

7 = <v(Jh> > = n kT P ' 

. 1 , > <V<7 S > 

2tt ; 2irkT * 



(8) 
(9) 



where k is the Boltzmann constant, p and T are the 
foreign-gas pressure and temperature, <...> denotes an 
average over both the gas and the beam velocity distribu- 
tions. The shift and broadening cross sections are given 
by (see, for example, Ref. |35j1: 



cr s = 2n I bsm((f>(b)) db, 
'o 



(10a) 



a b = 2n 6(1 - cos(0(&)) db. (10b) 
Jo 

Here, <j)(b) is the induced phase shift which depends on 
the form of the interaction potential. These potentials 
are usually expressed as Vfc=C^/r™ for the k th excited 
state with the van der Waals (n=6) being the primary 
long range interaction for neutral atoms. For this poten- 
tial, the collisional cross sections for a transition between 
arbitrary states a and b are given by 



<<7 6 >= 2tt(0.602) 



MCI 1x2/5 



<v> 



, (Ha) 



«j s >= sign(^) 2^(0.438) (j^j ^ , (Hb) 

where C' 6 = (Cg - C|)//L 

As illustrated by this example, a s /2a b has a value of 
0.364, independent of C 6 . In fact, within the adiabatic 
impact approximation, a unique ratio is predicted for 
each interaction potential which can be used to deter- 
mine the dominant form of the interaction. Since cross 
sections can have a dependence on the relative velocities, 
the quantities <<rv> for both broadening and shifts ex- 
tracted from our experimental results are tabulated in 



Table E] as well as their ratio j^. The sign convention 
for the shift cross sections has been chosen such that a 
positive frequency shift for the 3.1-MHz transition (where 
level B has higher energy than level A) gives a positive 
cross section. The measured ratios vary widely for differ- 
ent gases. In foreign gas perturbations, this is a common 
situation since in addition to long range attractive inter- 
actions, short range repulsive forces can be important. 
This is especially true for H2, He and Ne which do not 
interact strongly via the van der Waals interaction |37| . 
Our results are in agreement with this trend: H2 and 
Ne cause no shifts, while He induces an opposite shift 
as a result of dominant repulsive forces. On the other 
hand, all three gases show broadening comparable to the 
other gases. This can be understood by observing the 
integrands of Eqs. (|10fl . While the integrand for the shift 
cross section oscillates around zero, the integrand for the 
width cross section is always positive, yielding a nonzero 
broadening even for short-range interactions. 



VI. IMPLICATIONS FOR THE a- VARIATION 
EXPERIMENT 



Collisional shifts are one of the most important sys- 
tematics for an experiment searching for the temporal 
variation of a using Dy since they mimic frequency shifts 
that would be induced by variation in a. The a depen- 
dence of different isotopes and hyperfine-level manifolds 
of odd isotpes are the same [38|. As a result, the en- 
tire manifold of level B shifts identically with respect to 
the manifold of level A. Monitoring different transitions 
where energy of level A is larger than the energy of B, 
and vice versa provides a powerful tool in detecting and 
eliminating certain systematics since the shift induced 
by a variation should have the same magnitude but op- 
posite sign for these transitions. However, as seen in 
Table this behavior is identical to that of the shifts 
induced by collisions, and therefore, gas pressures must 
be kept stable to 1 part in 10 4 for a background pressure 
of ~ 1 /iTorr, or lowered to < 3 x 10~ 10 Torr to achieve 
the desired mHz-level stability over a year. 

This conclusion is in agreement with the estimates pro- 
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vided in Ref. ^l|. However, it was also argued in Ref. |Tl| 
that a suppression of collisional shifts might exist due to 
the fact that the transitions are nominally between f and 
d inner-shell electrons and that inelastic quenching might 
remove atoms experiencing collisions from levels A or B 
so that they would not contribute to the signal. Our 
present results do not reveal such suppression of colli- 
sional shifts. 

With the current experimental setup, the operating 
pressure is ~ 5 //Torr, mostly dominated by hydrogen 
and nitrogen. A new vacuum chamber that can achieve 
1CP 10 Torr with a Dy oven operating at 1500 K is under 
construction. In the mean time, data that is being ac- 
quired with the current chamber will be corrected using 
the measured collisional shift coefficients to achieve an es- 
timated stability of 2 Hz over a year, which corresponds 
to sensitivity of \a/a\ ~ 10~ 15 /yr. This sensitivity is 
comparable to the best current laboratory limits pol |22|| . 

VII. CONCLUSION 

In summary, we have reported collisional shift and 
broadening rates, and cross sections for radio-frequency 
transitions between excited states of 162 Dy and 163 Dy in 
the presence of 10 to 80 /xTorr of H 2 , N 2 , 2 , He, Ne, Ar, 



Kr, and Xe. In addition, we have placed upper bounds on 
quenching cross sections using the pressure dependence 
of the signal amplitude. A frequency-modulation tech- 
nique used for the search of the temporal variation of 
the fine-structure constant was utilized to measure the 
small perturbations induced at such low pressures. The 
extracted cross sections are similar to cross sections for 
optical transitions 

To improve our understanding of the interaction po- 
tentials, several avenues may be pursued in the future. 
As seen in Eqs. (fTTjl . the possible velocity dependence of 
cross sections can be exploited by varying the gas tem- 
perature. Moreover, a careful study of the line shape far 
off resonance would be informative due to strong temper- 
ature dependences and satellite-line formations at higher 
pressures |37| . 
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